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UNIVERSAL  RATE  MODEL  SELECTOR:  A  METHOD 
TO  QUICKLY  FIND  THE  BEST-FIT  KINETIC  RATE  MODEL 
FOR  AN  EXPERIMENTAL  RATE  PROFILE 


1.  INTRODUCTION 

A  kinetic  experiment  (e.g.,  time  vs  intensity  or  concentration  measurements) 
requires  a  kinetic  rate  model  to  fit  the  experimental  data  points.  In  the  literature,  there  are  many 
kinetic  rate  models  to  choose  from  ( 1 ,  2),  including  the  fractional  power  (3),  pseudo-first-order 
(4),  pseudo-second-order  (5),  Elovich  (6),  Avrami  (7),  and  intraparticle  diffusion  (6)  rate 
equations  to  name  a  few. 

A  single  fitting  algorithm  (kinetic  rate  model)  for  a  reaction  does  not  usually 
address  multiple  regions  of  an  experimental  rate  profile,  such  as  the  initial  or  end  (saturation) 
stages  where  deviations  or  discrepancies  may  be  found.  Experimental  data  versus  model-fitting 
discrepancies  at  the  initial  stage  of  a  rate  profile  analysis  are  common  because  the  independent 
y  values  are  generally  located  on  the  steep  portion  of  the  kinetic  rate  decay  curve.  This  is  also  a 
concern  for  the  later-stage  data  points  because  a  greater  degree  of  scatter  is  common.  Different 
chemical  and  physical  mechanisms  may  occur  during  an  experiment  to  cause  these  differences. 
The  fitting  algorithm  also  provides  the  best-fit  parameters  (degree  of  correlation)  to  produce  a 
minimum  error-squared  value.  As  the  number  of  fitting  parameters  in  a  rate  model  increases,  a 
higher  degree  of  correlation  is  usually  observed.  However,  over-fitting  may  occur,  and  the  fitting 
parameter(s)  may  lose  meaning  or  association. 

In  general,  a  kinetic  rate  model  can  be  divided  and  classified  according  to  the 
number  of  rate  constants.  It  should  be  noted  that  the  majority  of  traditional  rate  models  need  to 
be  converted  into  time  versus  relative  intensity  (or  relative  concentration  or  percent  conversion) 
so  the  values  can  be  directly  compared. 

The  concept  of  forcing  or  fitting  a  particular  kinetic  rate  model  onto  a  complete 
data  set  may  overlook  salient  or  significant  alternate  trends  in  a  portion  of  the  complete  record  of 
data.  Herein,  we  constructed  a  blueprint  (platform)  set  of  graphs  that  contains  eight  traditional, 
widely  used  kinetic  rate  model  curves  as  the  platform  for  the  Universal  Rate  Model  Selector 
(URMS).  The  URMS  was  tested  using  surface-enhanced  Raman  spectroscopy  (SERS)  of 
thiophenol  adsorption  onto  a  gold  surface.  SERS  experimental  data  sets  of  thiophenol,  which 
consisted  of  different  experimental  conditions  and  parameters  ( 8 ),  were  overlaid  directly  onto  the 
URMS  blueprint  platform  (eight  kinetic  rate  curves).  A  visual  observation  showed  where  the 
normalized  experimental  data  points  most  closely  associated  with  a  particular  rate  curve.  No 
fitting  or  calculations  were  performed  in  the  fit  between  experimental  data  and  the  URMS. 
Instead,  a  visual  analysis  was  conducted.  High  degrees  of  correlation  were  observed  with  the 
normalized  experimental  data  points  and  selected  URMS  curves.  Therefore,  the  experimental 
data  points  were  used  to  select  the  rate  model  curves  that  showed  the  closest  association. 
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2. 


MATERIALS  AND  METHODS 


Experimental  thiophenol  rate  profiles  (time  vs  SERS  signal)  at  four  different  pH 
values  and  six  different  temperatures  (total  of  25  rate  profiles)  were  acquired  using  commercial 
nanostructured  gold  SERS  substrates  (8).  Thiophenol  is  a  model  system  for  adsorption  rate 
studies  of  molecules  on  noble  metal  surfaces  due  to  the  strong  affinity  of  the  -SH  group  for  these 
metals.  Figure  1  shows  all  25  experimental  thiophenol  rate  profiles  at  four  different  pH  values 
(four  panels  in  Figure  1)  and  at  different  temperatures.  There  are  six  rate  profiles  at  pH  2  and 
temperatures  of  1.7,  13.7,  25.2,  34,  42.5,  and  54.1  °C.  There  are  six  rate  profiles  at  pH  4  at 
temperatures  of  0.7,  3.7,  13.3,  26.4,  42.6,  and  54.5  °C.  There  are  six  rate  profiles  at  pH  6  at 
temperatures  of  1.5,  12.8,  25.7,  34.2,  42.5,  and  55.3  °C.  There  are  seven  rate  profiles  at  pH  10  at 
temperatures  of  1.7,  3.4,  12.5,  25.2,  32.6,  42.8,  and  53.7  °C.  Different  rate  profiles  have  different 
maximum  intensities,  which  may  be  attributed  to  the  variations  on  the  SERS  substrate,  Klarite, 
slides  (Renishaw  Diagnostics  Limited,  Glasgow,  U.K.).  Intensity  spikes  may  be  due  to 
(1)  variations  between  Klarite  slides,  (2)  system  performance  over  a  period  of  time,  (3)  pH 
control,  (4)  temperature  control,  and  (5)  environmental  noise. 
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Figure  1.  Twenty-five  experimental  thiophenol  rate  profiles  at  pH  2,  4,  6,  and  10  and  at  various 

temperatures. 
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3. 


RESULTS  AND  DISCUSSION 


3.1  Calculation  of  Effective  Intensity 

The  first  step  was  to  scale  and  normalize  the  SERS  intensity  values  and  convert 
them  to  experimental  relative  intensities  (Rex P  =  0-1)  using  the  experimental  minimum  and 
maximum  intensity,  /min  and  /max,  respectively,  as  follows: 

Rex  p  —  (/  /min)/ (/max  —  /min)  ( 1 ) 


where  ReXp  is  the  experimental  relative  intensity,  based  on  the  experimental  maximum  intensity. 

The  question  then  arose  as  to  how  to  scale  the  intensity  values  (/)  and  determine 
the  proper  /min  and  /max  values.  For  example,  the  /max  value  or  the  average  maximum  intensity 
value  at  the  saturation  or  plateau  region  could  be  used  (Figure  2).  Equation  1  was  modified  as 
follows: 


Rp  —  Radj  X  Rexp  (2) 

where  Rp  is  the  predicted  relative  intensity,  based  on  the  predicted  maximum  intensity,  /m axp,  and 
Radj  is  the  adjusted  /max  intensity  factor 


R'dd]  —  (/max  —  /min)/(/maxp  —  /min)  (3) 

=  (/max  —  /min)/ (/max  +  A  /  —  /min) 

and 

/m  axp  =  /max  +  A /  (4) 

Radj  is  necessary  because  this  term  takes  into  consideration  whether  an  experiment  was 
performed  to  completion  (into  the  plateau  or  saturation  phase,  e.g.,  Figure  1  at  pH  2,  42.5  °C  and 
at  pH  6,  1 .5  °C)  or  whether  the  experiment  was  not  allowed  to  proceed  to  completion.  If  an 
experiment  is  stopped  before  it  reaches  the  saturation  or  plateau  phase,  then  that  /max  (apparent 
maximum  intensity)  is  actually  lower  than  the  true  maximum  (predicted  maximum  intensity). 
Therefore,  the  A /  term  is  the  difference  between  /max P  and  /max.  Examples  of  this  situation  are 
shown  in  Figure  1  at  pH  2,  1.7  °C  and  at  pH  6,  25.7  °C.  Therefore,  Rp  at  time  it)  is  the  adjusted 
Radj  intensity  multiplied  by  ReXp. 
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Figure  2.  Experimental  rate  profiles  at  pH  6.  The  /max  and 
average  maximum  intensity  values  are  indicated  with  arrows. 


Figure  3  shows  all  25  experimental  thiophenol  rate  profiles  at  pH  2,  4,  6,  and  10 
and  at  different  temperatures  after  the  intensity  values  were  converted  to  effective  intensities  ( Rp ) 
using  eq  2.  The  rate  profiles  at  pH  10  (Figure  3,  bottom  panel)  show  more  temperature 
dependency  than  the  rate  profiles  at  the  other  pH  values.  Rate  profiles  at  pH  2  show  the  least 
temperature  dependency,  but  the  profile  at  13.7  °C  is  questionable.  Rate  profiles  at  pH  6  show 
strong  temperature  dependency  at  temperatures  12.8  and  1.5  °C,  and  rate  profiles  at  pH  4  show 
some  degree  of  temperature  dependency. 
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Figure  3.  Twenty-five  experimental  thiophenol  rate  profiles  at  pH  2,  4,  6,  and  10  and  at  various 
temperatures  after  the  intensity  values  were  converted  to  effective  intensities  using  eq  2. 
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3.2 


Kinetic  Rate  Models 


Metal-organic  complex  mechanisms  of  adsorption  often  involve  a  chemical 
reaction  between  functional  groups  on  a  sorbent  surface  with  the  metal  ions,  thus  forming  metal- 
organic  complexes.  This  produces  a  cation-exchange  reaction  due  to  the  high  cation-exchange 
capacity  of  the  sorbent.  In  addition,  other  mechanisms  may  be  involved  including  mass-transport 
processes  such  as  transport  in  the  bulk  of  the  liquid  phase,  diffusion  across  the  liquid  film  near 
the  surface,  and  diffusion  in  macro-  or  micropores.  Adsorption  kinetics  have  been  extensively 
studied.  Usually,  the  adsorption  rate  is  very  rapid  at  the  beginning  of  the  process  and  then 
becomes  slower  as  equilibrium  is  approached.  Such  kinetics  are  often  described  by  a  first-order 
kinetic  reaction  (9-11)  and  are  sometimes  interpreted  as  a  combination  of  two  or  three  different 
mechanisms  (12-18).  McKay  and  Allen  (19,  20)  studied  the  sorption  of  dyes  onto  peat  and 
developed  mass-transport  diffusion  models  that  were  used  successfully  to  predict  the  dye 
concentration  versus  time-decay  curves. 

Table  1  shows  some  of  the  common  kinetic  rate  models  with  one  rate  constant 
that  are  found  in  the  literature  (7,  21,  22).  Table  2  shows  some  of  the  common  kinetic  rate 
models  with  two  rate  constants  that  are  found  in  literature  (1-7,  21,  22). 


Table  1.  Kinetic  Rate  Models  with  One  Rate  Constant 


No. 

Name 

Rate  Equation 

1st  tso 

2nd  /so 

1 

First-order  model 

•  Film-diffusion 
model, 

•  Mass  transfer 
rate  or  linear 
driving  force 
(LDF),  and 

•  Pseudo-first- 
order  model 

ln(l  -  Rp)  =  ~(knt) 

ln(l  -  Rp)  =  -(at) 

Ei/RT  = 

lnfeo)  -  ln[-ln(0.5)] 

ln(l  -  Rp)  = 
(r/t5o)ln(0.5) 

2 

Second-order  model 

Pseudo-second- 
order  model 

[1/(1  -RP)]  = 

1  +  (fc^)/max 

EORT  = 

lnfeo)  -  ln[(0.5_1  -  1)] 

(1-^)-*  = 

1  +  (t/t50)  [0.5-1  -  1] 

3 

nth-order  model 

(l-/?p)(1^  = 

1  +  ( kn  t) 

(jl  —  1)  /max 

Ei/RT  = 
ln(t5o)  - 

ln[(0.5(1-n)  -  l)/(n  -  1)] 

(1  -  Rp)(1_")  = 

1  +  (t/tso)  [0.5(1_n)  -  1] 

4 

Film  diffusion 
model,  first-order 
model  (Dumwald- 
Wagner) 

ln(l  -RP2)  = 

-(kit) 

Ei/RT  =  ln(t5o) 

-  ln[-ln(l  -  0.52)] 

ln(l  -RP2)  = 

(t/t5o)ln(l  -  0.52) 

(table  continued) 
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Table  1.  Kinetic  Rate  Models  with  One  Rate  Constant  (Continued) 


No. 

Name 

Rate  Equation 

1st  *50 

2nd  tso 

5 

First-order 
autocatalytic  reaction 
model 

Rp  +  a  = 

[1  +  a]/ 

[1  +  (l/a)e-(feXP)] 

Ei/RT  =  ln(t5o)  - 
ln[-ln(0.5a)/(0.5  +  a)] 

ln[a(l  -  Rp)/(Rp+a)]  = 
(t/f5o)ln[a(l  -  0.5)/ 

(0.5  +  a)] 

6 

Avrami  rate  model 

ln(  1  -  Rp)  = 

-m 

Ei/RT  = 

nln(f5o)  -  ln[-ln(0.5)] 

ln(l  -  Rp)  = 

(t/tsoT  ln(0.5) 

Nonconverging  Kinetic  Rate  Models  (R  oo  when  t  oo) 

7 

Fractional  power 
model 

RP  =  ( ktb ) 

Ei/RT  = 

Mn(f5o)  -  ln(0.5) 

Rp/0.5  =  (r/tsof 

8 

Intraparticle  diffusion 
model  (Weber- 
Morris) 

Rp  X  /max  = 

Cd  +  kd  t0  5 

Ei/RT  = 

0.51n(t5o)  - 
ln(0. 5/max  —  Cd) 

(Rp  —  Cd//max)  = 

(r/r50)0  5  (0.5  -  Cd//max) 

9 

Bioadsorption: 

Elovich  model 

Rp  =  ln(l  + 

debei)/ (bel max) 

Ei/RT  = 

ln(f5o)  -  ln[(ea5W 
-  1  )/be)] 

When  abtso  »  1 

Rp  -  0.5  = 

ln(//f50  )/ (bel max) 

where 

/50  is  time  when  the  reaction  reaches  a  relative  intensity  of  0.5;  k  is  rate  constant;  t  is  time;  R  is  ideal  gas  constant;  T  is 
temperature;  a  is  mass  rate  diffusion  rate  or  LDF  effective  mass  transfer  coefficient;  E\  is  activation  energy  associated 
with  rate  constant  ki;  kf  is  film  diffusion  constant;  a  is  catalyst  ratio  equal  to  /o//max;  P  =  (1  +  a)  Imax  ;  n  is  a  constant;  b  is 
a  constant  with  b  <  1 ;  kdis  rate  of  intraparticle  diffusion-controlled  adsorption  constant;  Cd  is  0  when  diffusion  is 
limiting  step;  ae  is  the  initial  adsorption  rate  (mg/g-h);  and  be  is  desorption  constant  (g/mg). 
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Table  2. 1 

ECinetic  Rate  Models  with  Two  Rate  Constants 

No. 

Name 

Rate  Equation 

2nd  tso 

1 

Homogeneous 
catalytic  reaction 
model 

ln(l  -  Rp)  =  ~{k\  +  kiSi)t 
ln(l  —  Rv)  =  —(kobst ) 

ln(l-Rp)  =  (r/r5o)ln(0.5) 

2 

Shifting  order 
model  1: 

Low  to  high 

— ln(  1  —  Rp)  +  fe/maxRp  = 

(kit) 

kll max  —  ln(  1  —  Rp)l Rp  = 

(asi  )(t/tso)/Rp 

3 

Shifting  order 
model  2: 

High  to  low 

ln(as2/[as2  -  fe7maxRp])  = 

(kit) 

ln(  1  —  fc/max^p/ (Xs2)  = 

(f/t5o)ln(l  —  0.5fo/max/ots2) 

4 

Double  exponential 
model  (DEM) 

( W  ilezak-Keinath) , 
two  adsorption  sites 
(slow,  rapid) 

d-Rv)  = 

aie-*T  +  (1  -  ai)e-k2( 

NA 

5 

Monomolecular 

consecutive  model 

OR 

In-series  reaction 
(n  =  2)  model 

(1-Rp)  = 
ae-*T  +  (1  -  a)e-*2f 

NA 

where 

S  is  available  sites  (S2  is  type  2  sites);  k0bs  is  observed  (effective)  rate  constant;  NA  is  not  applicable. 
kobs  =  k\  +  kiSi ;  asi  =  0.5 fe  -  ln(0.5)//max ;  aS2  =  k\  +  fe/max;  a\  =  fraction  of  site;  a  =  fa/(fa  -  k\) 

3.3  Universal  Kinetic  Rate  Platform  Development 

Kinetic  rate  models  range  from  pure  chemical  reactions  to  mass  transfer 
(diffusion)  (Tables  1  and  2).  The  rate  equations  in  Tables  1  and  2  were  converted  to  be 
independent  of  the  time,  intensity,  temperature,  and  activation  energy  parameters.  This  new 
approach  offers  a  unique  single  curve  for  each  rate  model  equation.  In  this  format,  there  is  no 
need  to  obtain  or  evaluate  kinetic  rate  constants  or  mechanisms.  The  experimental  data  can  be 
mapped  visually  onto  the  transformed  rate  models  that  provide  a  best  fit.  The  fit  between  the  data 
points  and  the  model  curves  requires  no  assumptions,  calculations,  or  prior  information  about  the 
experimental  conditions. 

All  rate  models  in  Tables  1  and  2  are  in  the  time  and  intensity  domains,  and  the 
time  parameter  needs  to  be  converted  into  dimensionless  units.  To  do  that,  we  used  the  tits o 
parameter.  The  tso  is  the  time  when  the  reaction  reaches  a  relative  intensity  of  0.5  (Rp  =  0.5),  or 
the  time  t  when  the  reaction  reaches  50%  conversion.  There  are  two  approaches  for  constructing 
the  tso,  as  shown  in  Table  1.  The  first  tso  approach  (1st  tso)  was  applied  to  the  rate  models  in 
Table  1,  and  the  results  are  shown  under  the  1st  tso  column  in  Table  1.  Table  1  shows  a  general 
linear  relationship  between  ln(tso)  and  1/7’,  where  the  slope  equals  E\/nR  and  n  =  0.5  for  rate 
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model  no.  8  (Weber-Morris),  n  is  constant  for  rate  model  no.  6  (Avrami  model),  and  n  =  1  for  all 
other  models. 


The  1st  tso  approach  is  a  convenient  model  used  to  produce  a  rough  estimate  on 
the  observed  E\/R  values  from  the  experimental  rate  profiles  at  different  temperatures.  However, 
the  1st  tso  approach  is  dependent  on  temperature  and  E\/R. 

The  second  tso  approach  (2nd  tso)  was  applied  to  the  rate  models  in  Table  1. 

Table  1  shows  that  (t/r so)"  is  a  function  (f)  of  Rp  or  f(Rp)  with  n  =  0.5  for  rate  model  no.  8;  n  is 
constant  for  rate  model  no.  6,  and  n  =  1  for  all  other  models.  The  2nd  tso  approach  appears  to  be 
superior  to  the  1st  tso  approach  because  the  2nd  tso  approach  converts  the  traditional  kinetic  rate 
models  so  that  they  are  independent  of  time,  intensity,  temperature,  rate  constant,  and  activation 
energy. 


Table  1  shows  that  each  kinetic  rate  model  has  a  single,  unique  curve  in  the  (t/tso) 
versus  Rv  plot.  Because  (t/tso)  =  1  is  an  inversion  point  in  the  kinetic  rate  models,  Rp  was  replaced 
with  i?p(l  -  Rp)  to  show  the  initial  stage  and  the  later  (close-to-completion)  stage  curvatures  for 
each  rate  model.  The  kinetic  rate  curves  generated  by  (t/tso)  on  the  x  axis  and  Rp(  1  -  Rp)  on  the 
y  axis  become  the  URMS  platform. 

Results  from  the  rate  models  in  Table  2  were  placed  into  the  2nd  tso  approach  as 

listed. 


Figure  4  shows  the  URMS  platform  where  the  kinetic  rate  models  have  no 
dependence  on  temperature,  time,  and  Ei/R.  The  four  different  temperatures  (1.7,  10,  25,  and 
43  °C)  produced  identical  curves  in  the  URMS  platform  for  each  respective  rate  model,  which 
showed  kinetic  rate  model  temperature  independence.  For  the  one  rate  constant  models 
(Table  2),  there  was  no  dependency  on  E\/R.  This  was  also  true  for  the  two  rate  constant  kinetic 
models  when  |A£i2l/R  was  kept  constant.  AEn  =  £2  -  £),  where  £2  and  E\  are  the  activation 
energies  associated  with  rate  constants  fe  and  k\,  respectively. 
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T(°C)  =  1.7  T(°C)  =  10 


T("C)  =  25  T{°C)  =  43 


Figure  4.  One  rate  constant  models  in  the  URMS  platform  that  are  independent  of  temperature, 

time,  and  Ei/R. 


Figures  5  and  6  show  the  effects  of  changing  the  Ei/R  and  \AEn\/R  values, 
respectively,  for  the  two  rate  constant  kinetic  models.  Figure  5  shows  identical  plots  for  the 
URMS  platform  when  \AEu\IR  values  are  kept  constant  (\AEu\/R  =  55)  while  changing  the  E\/R 
values.  But  when  \AEn\IR  increases  (Figure  6),  the  curve  broadens  significantly  on  the  right  side 
of  the  plot  (t/tso  =  1),  regardless  of  whether  E\/R  is  held  constant  or  not.  Each  \AEn\/R  value 
produces  a  unique  curve  in  the  URMS  platform.  The  percentage  of  (\AEn\/Ei)  appears  to 
strongly  affect  the  later  stage  degree  of  widening,  especially  above  50%. 

Figure  6  shows  a  general  trend  on  URMS  curves  containing  two  rate  constants. 
For  each  \AEn\/Ei  value,  there  exists  a  unique  URMS  curve  with  the  two  rate  constant  kinetic 
rate  model. 
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The  dotted  and  dashed  curves  are  reference  curves  with  fixed  Ei/R  values  of  3870.  It  is  clear  that 
the  blue  and  green  curves  did  not  change  as  compared  with  the  reference  curves  when  \AEn\/R 
values  were  kept  constant  (55),  while  the  Ei/R  values  were  changed.  Cons.,  consecutive  (in¬ 
series)  rate  model. 
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URMS  platform.  When  \AEn\/R  was  increased,  the  DEM  (a)  showed  more  significant  widening 
than  the  Cons,  (b)  for  the  same  (\AEn\/Ei)  percentage. 


Figure  7  shows  the  URMS  curves  with  nth-order  rate  models  at  n  =  0.25,  1,  and  2; 
homogeneous  catalytic  rate  model  (pseudo-first-order);  DEM  at  two  different  \AEn\  values 
(55  and  155);  and  consecutive  rate  models  at  two  different  \AEn\  values  (55  and  155).  We  chose 
these  rate  models  because  they  represent  the  commonly  used  traditional  kinetic  rate  models 
found  in  the  literature  (3-11,  21,  22).  Rp  =  0-1  when  the  time  (t)  goes  from  0  to  oo. 
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Figure  7.  URMS  models  at  n  =  0.25,  1,  and  2.  Also  shown  are  the  DEM  and  Cons,  at  selected 

AE\  values. 


We  tested  the  URMS  platform  (Figure  7)  with  experimental  SERS  thiophenol  rate 
profiles  at  pH  10  (Figure  8).  The  six  experimental  rate  profiles  at  1.7,  12.5,  25.2,  32.6,  42.8,  and 
53.7  °C  were  evaluated  visually  to  see  which  rate  model  had  the  highest  degree  of  correlation 
with  the  experimental  data  points.  The  URMS  platform  offers  a  unique  type  of  evaluation,  which 
is  a  comparison  of  each  experimental  rate  profile  with  any  changes  with  varying  temperatures. 
Figure  8  shows  that  a  first-order  rate  model  is  generally  a  best  fit  for  the  experimental  points. 
This  behavior  was  also  noted  by  Biggs  and  colleagues  (23),  where  the  SERS  peak  growth 
appeared  to  follow  the  time-dependent  Langmuir  kinetics  (first-order  rate  model). 

There  were  some  deviations  from  a  first-order  rate  model  at  the  later  stage  when 
t/tso  >  2.5.  Was  this  because  (1)  the  Rp  values  were  very  similar  to  Rexp  or  (2)  there  was  a  change 
in  the  mechanism  from  a  first-  to  second-order  reaction  rate  model?  Biggs  and  colleagues  (23) 
noted  that  the  DEM  behavior  was  consistent  with  a  gradual  reordering  of  the  adsorbed  layer  as 
the  surface  approached  saturation  coverage.  This  is  a  phenomenon  previously  observed  in 
benzylthiol  and  other  organothiol  adsorption  studies  (24-26). 

Also,  there  was  a  change  in  the  initial  stage  of  the  profiles  when  t/tso  <  1 .0.  Some 
of  the  experimental  data  points  were  closer  to  a  second-order  rate  model,  whereas  other  data 
points  were  closer  to  a  first-order  model. 
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a  first-order  or  homogeneous  catalytic  reaction. 


We  also  tested  the  URMS  platform  in  Figure  7  with  a  second  set  of  experimental 
SERS  thiophenol  rate  profiles  at  pH  6  (Figure  9).  The  six  experimental  rate  profiles  at  1.5,  12.8, 
25.7,  34.2,  42.5,  and  55.3  °C  were  evaluated  visually  to  see  which  rate  model  had  the  highest 
degree  of  correlation  with  the  experimental  data  points.  Figure  9  shows  that  no  rate  model 
adequately  fit  the  experimental  data  points.  The  consecutive  rate  model  showed  a  low 
correlation. 


0.25 

0.2 

S' 

±0.15 

oc 

0.1 

0.05 

0 


pH  =  6 


(*/  *5o) 


DEM 


\AE\=  55&155 


lst-order  =  Homogenous  Catalytic  Rx 


Figure  9.  No  single  rate  model  appears  to  fit  the  experimental  data. 


We  tested  the  URMS  platform  in  Figure  7  with  a  third  set  of  experimental  SERS 
thiophenol  rate  profiles  at  pH  4  (Figure  10).  Figure  10  shows  that  the  six  experimental  rate 
profiles  at  -0.7,  3.7,  13.3,  26.4,  42.6,  and  54.5  °C  were  evaluated  visually  for  the  rate  model  that 
had  the  highest  degree  of  correlation  with  the  experimental  data  points.  No  single  rate  model 
provides  an  adequate  fit  to  the  experimental  data  points.  The  consecutive  rate  model  had  a 
substantially  low  degree  of  correlation  with  the  data  points.  It  is  interesting  that  Figures  9  and  10 
showed  similar  experimental  data  point  profiles.  This  was  significant  because  at  pH  10,  the  rate 
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profiles  were  different  from  those  at  pH  6  and  below.  This  behavior  was  also  explained  by  Biggs 
and  colleagues  (23),  where  the  growth  in  intensity  of  the  422  cm-1  peak  was  best  described  by 
the  DEM  two  rate  constant  model.  The  422  cm-1  curve  deviation  from  first-order  Langmuir 
kinetics  was  likely  related  to  its  strong  dependence  on  the  azz  tensor  value,  that  is,  its 
electromagnetic  enhancement  was  very  sensitive  to  surface-layer  orientation  (27). 


OC 

A  0.15 


pH  =  4 


(t/t50) 


Ist-order  =  Homogenous  Catalytic  Rx 


Figure  10.  A  comparison  of  the  rate  models  with  the  experimental  data  at  pH  4. 


We  tested  the  URMS  platform  in  Figure  7  with  a  fourth  set  of  experimental  SERS 
thiophenol  rate  profiles  at  pH  2  (Figure  11).  The  six  experimental  rate  profiles  at  1.7,  13.7,  25.2, 
34,  42.5,  and  54.1  °C  were  evaluated  visually  to  see  which  rate  model  had  the  highest  degree  of 
correlation  with  the  experimental  data  points.  Figure  1 1  shows  that  no  combination  of  rate 
models  adequately  fits  the  experimental  data  points. 


Figure  1 1 .  A  comparison  of  the  rate  models  with  the  experimental  data  at  pH  2. 
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It  is  significant  to  note  that  from  pH  10  to  2  (Figures  8-11),  the  experimental 
profiles  showed  a  tendency  to  become  narrower.  Furthermore,  the  initial  stage  of  the 
experimental  profiles  (0  <  t/tso  <  1)  showed  a  delay  in  the  initial  rise  as  the  pH  value  decreased. 


4.  CONCLUSIONS 

The  URMS  presents  a  fast  and  easy  way  for  a  visual  analysis  of  complex 
experimental  rate  profiles  that  consist  of  many  data  points  obtained  at  different  operating 
conditions  (temperature  and  pH  values).  The  URMS  platform  helps  to  visually  evaluate  the 
experimental  rate  profiles  to  determine  the  best  kinetic  model  that  exhibits  a  relatively  high 
degree  of  correlation.  Any  changes  that  may  occur  during  the  initial  and  later  stages  of 
experimentation  (saturation),  as  well  as  how  the  rate  profile  changes  with  operating  conditions 
such  as  temperature  and  pH,  are  easily  observed  with  the  URMS.  The  kinetic  rate  model 
equations  were  solved  without  considering  the  major  parameters  such  as  temperature,  rate 
constant,  activation  energy,  time,  and  intensity. 

Future  efforts  will  address  rate  models  with  equal  to  or  greater  than  three  rate 
constants  (three  activation  energies;  £i,  £2,  and  £3).  URMS  models  will  be  devised  for  each  pair 
of  \AEn\/Ei,  |A£i3l/£i,  and  \AE23VE2  values.  A£i3  =  £3  -  E\  and  AE23  =  £3  -  £2,  where  £3  is  the 
activation  energy  associated  with  rate  constant  ks. 
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